Lipid membrane instability and poration driven by capacitive charging 
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A new model for the interaction of an electric pulse with a lipid membrane is proposed. Using this 
model we show that when a DC electric pulse is applied to an insulating lipid membrane separating 
fluids with different conductivities, the capacitive charging current through the membrane drives 
electrohydrodynamic flow that destabilizes the membrane. The instability is transient and decays 
as the membrane charges. The bulk conductivity mismatch plays an essential role in this instability 
because it results in a different rate of charge accumulation on the membrane's physical surfaces. 
Shearing stresses created by the electric field acting on its own induced free charge are non-zero as 
long as the charge imbalance exists. Accordingly, the most unstable mode is related to the ratio of 
membrane charging time and the electrohydrodynamic time. 

PACS numbers: 47.20.Ma, 47.57.jd, 87.16.dj 



S 1 
6 

4^ 

O 

Oh 



> 

o 

q 

in 

o 
o 



X 



Cells and cellular organelles are enveloped by lipid 
membranes, whose integrity is crucial for the cell viabil- 
ity. External electric fields can cause membrane break- 
down and cell death. However, a controlled application of 
electric pulses can induce transient pores in the cell mem- 
brane, which reseal after the pulse is turned off. This re- 
versible electroporation is of great interest for biomedical 
technologies such as gene transfection because the pores 
enable the delivery of exogenous molecules (drugs, pro- 
teins, DNA) into living cells pQ. However, controlled elec- 
troporation remains elusive because the physical mecha- 
nisms underlying membrane response to electric fields are 
poorly understood [2HI]. Here we propose a new model 
for the interaction of an electric pulse with a lipid mem- 
brane. The model sheds new light on the electroporation 
process and is applicable to a broad range of biological 
interfaces. 

The lipid bilayer membrane is impermeable to ions 
and acts as a capacitor when an electric field is applied. 
Charges carried by conduction accumulate at the phys- 
ical surfaces and a potential difference across the mem- 
brane builds up. In the case of a non-conducting mem- 
brane subjected to a uniform DC electric field, the trans- 
membrane potential increases as 



Vrr, 
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(1) 



where t m is a characteristic charging time for the capac- 
itor and V, in the case of a planar membrane, is the 
applied potential difference. It is generally accepted that 
the membrane porates when the the transmembrane po- 
tential exceeds a critical value of about one Volt [5]. In a 
typical electroporation experiment [5J, t m ~ 1ms. This 
is comparable and in certain cases even longer than the 
duration of the electric pulse, yet poration may occur 
at voltages lower than the critical one. This observation 
has been attributed to the presence of an initial tension in 
the membrane. A possible alternate explanation is that 



the membrane becomes unstable even before the criti- 
cal voltage is reached. In this Letter we propose a new 
mechanism for such an instability. We examine a non- 
conducting membrane and show that upon application of 
a uniform DC electric field a transient induced charge ap- 
pears on the membrane. The electric field acting on this 
charge gives rise to membrane in-plane electric stresses, 
which move the lipids and adjacent fluids. The resulting 
electrohydrodynamic flow can enhance a small perturba- 
tion of the membrane shape. This mechanism is similar 
to the destabilization of a fluid interface due to a finite 
time for charge relaxation [7 . 

Although the theoretical modeling of the dynamics of 
lipid membranes in steady electric fields has received 
considerable attention [SHllj. we are not aware of any 
study that has examined the effect of the transient elec- 
tric field on the stability of the membrane. Our model is 
fundamentally different from earlier work, in which vari- 
ous mechanisms for electroporation have been examined: 
pore nucleation [T21H3], membrane thinning FHrll6| , neg- 
ative tension [H] or induced charge electro-osmotic flow 
QTJ [T71 HH]. Moreover, unlike recent electromechanical 
models [3 fTMTB] . which have considered a conducting 
membrane, we analyze an insulating membrane. 

Let us consider a planar non-conducting membrane 
formed by a charge-free lipid bilayer with dielectric con- 
stant e mm . The bilayer thickness is about d ~ 5nra, 
thus on macroscopic length scales the membrane can be 
regarded as a two-dimensional surface with capacitance 
C m = fimm/rf- The membrane separates two fluids: a su- 
perphase of viscosity \i\ , conductivity <j\ , and dielectric 
constant t\, and a subphase characterized respectively by 
fJ-2, &2i and £2- The membrane is subjected to a perpen- 
dicular electric pulse with magnitude E — V/2L, where 
V is the applied potential, and 2L is the distance between 
the electrodes. The problem is sketched in Figure [T] 

Upon application of the electric field, bulk phases be- 
come electroneutral on a very fast time scale given by 



t c ,k = £k/&k (where fc = 1, 2 denotes the top or bottom 
fluids). Accordingly, the electric potential is a solution 
of Laplace's equation and the equations of bulk fluid mo- 
tion have no electric terms, i.e. the electromechanical 
coupling occurs only at the interface. This is in essence 
the leaky dielectric model developed by G. I. Taylor [!§] . 
Note that even though electrokinetic transport is not ex- 
plicitly included, the effect of space charge is accounted 
for in an aggregate sense [2"0Il2"T] . 




FIG. 1: A sketch of the problem. 

The electrostatic problem is specified by V • E& =0 
and V x Efc = 0, which is equivalent to E^ = — W<f>k and 
V 2 0fe = with boundary conditions at the electrodes 
<t>\{L) = -V/2, fa(—L) = V/2. At the membrane inter- 
face, z = h(x, t), conservation of normal currents requires 

dVrr 
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dt 



=<7in • Ei = o^n • E2 , 



(2) 



where Ohmic conduction is assumed for the bulk fluids. 
V m = 02 — 01 is the transmembrane potential and n is the 
unit normal vector. At t = when the pulse is applied 
the potential is continuous, V m (t = 0) = 0. The differ- 
ence in bulk fluid conductivities gives rise to a transient 
imbalance in the charge densities on the two sides of the 
membrane (since charges are brought at a different rate) 



= n (eiV0i -£2V0 2 )- 



(3) 



Q is not to be mistaken with the charge of the capacitor; 
if the capacitor is fully charged Q = 0. 

The electric field stresses the interface, which responds 
by deformation. The fluid flow accompanying the inter- 
face motion is described by the Stokes equations, because 
at the length scales of the membrane undulations viscos- 
ity dampens fluid acceleration. The velocity, v^, and the 
pressure, pk, fields satisfy 



Vp fc =MfcV 2 v fe , V-v fe = 0. 



(4) 



At the electrodes the fluid satisfies the no slip bound- 
ary condition and the velocity is continuous across the 
interface. The interface moves with the velocity of the 
surrounding fluid, 

dh /•„.„, 2N 1 / 2 



The surface flow is area-incompressible, V s • v s =0. 

The coupling of the electric and hydrodynamics fields 
at the interface is reflected by the stress balance 



(p 2 - Pl )n + n-[[T hd + T el }} =f 



(6) 



where [[...]] denotes a difference between fluid 1 and 



fluid 2. 



The Maxwell stress tensor is T„ e ; — 



dt 



n 1 



d x h 2 ) 



vi = v 2 



(5) 



e (EiEj - l/2E 2 5ij), T^ d = /1 (piV 3 + djVi) is the viscous 
stress tensor, and f m are the membrane tractions. In the 
simplest model (Hclfrich-Canham) 

f m = [-K (2H 3 - 2HK + V 2 H) + EH] n + V S S , (7) 

where k is the bending rigidity, H and K are the mean 
and Gaussian curvatures. Here £ is the membrane ten- 
sion, a Lagrange multiplier which enforces the surface- 
incompressibility [2"2] . 

Next we proceed to analyze the linear stability of a 
planar membrane subjected to an electric pulse. The 
evolution of a membrane perturbation in the x-direction 
with a wavenumber q depends on the interplay of several 
physical process: the decrease of the effective interfacial 
charge (which occurs on a characteristic time scale cor- 
responding to the charging of the membrane capacitor, 
t m ), shape distortion (which occurs on electrohydrody- 
namic time scale, t e hd), and curvature relaxation on a 
time scale t K 

t m = C m L/a 2 , tehd = ^2/(e2-Eo) j *« = M2/(k<7 3 )- 

(8) 
In a typical electroporation experiment, the fluids are 
aqueous salt solutions, a ~ 10 _4 5/m, C m ~ 0.01F/m 2 , 
and Eq ~ IkV/cm. Hence, t m ~ t e hd ~ lms while 
t c ~ 1/is, i.e. both the electrohydrodynamic and the 
capacitor charging times are much longer that the bulk 
charge relaxation time, which justifies the quasi-static 
approximation. 

The base state of the problem is a flat interface, h = 0, 
stressed by a perpendicular electric field; any variable 
g is expanded in a series of the form g = g(°'(z,t) + 
g^\z,x,t) + ..., where g^ = J^ q 9q{ z , *) exp(igir). The 
problem solution involves the following steps. First, we 
solve the electrostatic problem to find the electric trac- 
tions exerted on the membrane. Second, we solve for 
the fluid flow needed to satisfy the stress balance, Eq. 
[6] Finally, the interface evolution is obtained from the 
kinematic condition, Eq. [5j For the membrane to de- 
form, the normal electric stress must be nonvanishing 
riJ } = E^E^K This requires that the base state elec- 
tric potential is non-uniform, which occurs only if there 
is current through the membrane; for a purely capaci- 
tive interface this is possible only during the period of 
charging. 

Henceforth, all variables are nondimensionalized. The 
characteristic time scale is the capacitor charging time 
t m , the length scale is L = V/Eq, and stresses are scaled 



by the viscous fluid stress \ii/t m . We define the conduc- 
tivity ratio R = cri/a^, permittivity ratio S = ei/e2, and 
viscosity ratio A = //i//«2- 

Base state: Flat membrane 

The solution for the transient electric field about a flat 
membrane yields 



■\ + {z 



b^ = \ + {z + l)RE^ 



l)E[°\t). 



(9) 



where e[ 0) (t) 



7(1 + .R) and a = R/(l+R). Hence, 



the effective surface charge and transmembrane potential 
are 



Q 



(0) 



(S-R)E[°\t), VW = l- e - 



(10) 



Initially the potential is continuous across the membrane, 
but the interface is charged (positively for R/S < 1, i.e., 
if the bottom fluid is more conducting, and negatively for 
R/S > 1). At times much longer than the capacitor time 
t m the membrane becomes fully charged, the electric field 
in the bulk vanishes and the transmembrane potential 
becomes equal to the applied potential. 

Fluctuating membrane: leading order analysis of the 
electric field effect on membrane undulations 

The solution of Laplace's equation for the electric po- 
tential in a bounded domain gives 



^ = ^(t)sinh[^-l)], 
<f>W=RA q (t)smh[q(z + l)]. 



(11) 



The coefficient A q (t) is determined from the conservation 
of current condition (Eq. [2]), which at this order couples 
to the amplitude of the interface perturbation h q (t) 



dA q qRcoth(q) 
~~dt~ + 1 + R 



A n 



l ~ R V*f) 



sinh(<7) (1 + R) dt 



(12) 

Note that the surface incompressibility requires v x {z = 
0) = 0. 

At this order in perturbation, the kinematic condition 
is dh/dt — v z (z — 0). Inserting the hydrodynamic and 
electric stresses in the normal stress balance (Eq. [6] and 
Eq. [7} yields 



1 + 2q 2 - cosh2g / , 3 , 

2 g + smh2 g (^ +Cg) ^( 13) 
-IL 2 (R 2 - S)A q (t)E[°\t) coshq) 



, . dh„ 

2(1 + A) ir 



where III = t m /t K (t K is defined based on L), II2 = 
tm/tehd) an d C = L 2 T, eq /K. We see that while the first 
term on the right hand side describes relaxation, the sec- 
ond term is destabilizing if E[ (t) =£ and there is a 
mismatch in the conductivities of the fluids. The tan- 
gential stress balance gives a nonuniform membrane ten- 
sion E,(i) = {R 2 - S)E[°\t)[A q (t)+E[°\t)h q (t)} which 
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FIG. 2: (Color online) Increasing shape undulation in response to 
a perturbation in the transmembrane potential. The arrows show 
the electric field at t = 0. The color coding reflects the magnitude 
of electric shearing stresses (red = zero shear stress). Interface 
profiles at t = 0, t = 0.5, and t = 2. Parameters as in Figure [3] 



influences the dynamics at next order. To illustrate this 
time dependent behavior, Eq. [12] and Eq. [13] are solved 
numerically for q — 4 and an initially flat interface (ini- 
tial condition small A4 7^ and hi = 0). In response to 
the electric pulse the membrane deforms; this is accom- 
panied with variations of the effective interfacial charge 
and shearing stress, see Figure [2] The tangential electric 
field is discontinuous (as a result of the presence of trans- 
membrane potential) producing vorticity and a behavior 
similar to the Kelvin-Helmholtz instability. The shearing 
stresses decrease in time; once they vanish, the growth 
of the interface perturbation stops and it begins to relax 
back to the flat, unperturbed state. Increasing the ap- 
plied voltage magnifies the amplitude of the perturbation 
in the shape of the membrane. 
Eq. 



[12] and Eq. 13 along with base-state electric 
field E[ '(t) describe the evolution of the membrane 
perturbation. This system can be written concisely as 
d t g = C(t)g where g = (h q (t),A q (t)) T . In general, a 
numerical solution of this system is required to investi- 
gate the stability of the membrane. Analytical results, 
however, can de derived if we adopt the frozen coefficient 
approximation. In this approach, the matrix C(t) is as- 
sumed constant allowing us to seek a solution of the form 
g(x, z, t) — g(z) exp[iga; + tut}. The growth rates Wi,2 are 
found as eigenvalues of C(t* ) for a fixed time t* . Although 
this approach is only an approximation, the stability re- 
sults are consistent with our numerical solution of the 
linear system of equations. 

In Figure [3] the dispersion relationship is shown for 
various frozen times at a fixed voltage of V = 0.14. As 
time progresses, the number of unstable wavenumbers 
decreases and eventually reaches zero. For small q the 
dispersion relationship seen in Figure [3] can be approxi- 
mated by 
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FIG. 3: Time evolution of the dispersion relation at fixed voltage. 
The dashed, dotted, and solid lines correspond to t* = 0, t* = 1, 
and t* = 2, respectively. The parameters are: R = 2, S = 1, E = 



lAkV/cm, k = 10" 19 J . 
<ti = 10 — 4 F/m.s, and C n 



= lO-Sj/m^ 
10- 2 F/m 2 . 



5.68 x 10" 



3 F/m, 



where 9{t*) = U 2 (E[ 0) (t*)) 2 'R, 1 /(l + A), K x = (R - 
l)(R 2 - 50/(1 + R) and <5(f ) = 9{t*)(l + R). The 
first term shows that the electric field acts as negative 
tension to destabilize the membrane; its magnitude di- 
minishes with time and eventually the restoring bending 
and membrane tension forces take control. Setting wi to 
zero in Eq. 14 allows us to identify a cutoff wave num- 



ber. The result confirms that instability can only occur if 
(R — l)(i? 2 — S) > 0. Even in the absence of permittivity 
mismatch (S = 1), the instability occurs if the bulk fluid 
conductivities differ. If S ^ 1, then (R - l)(R - y/S) 
must be positive to produce an instability for a nonzero 
electric field. 

We have also analyzed the electrohydrodynamics of a 
membrane modeled as two coupled monolayers [531 121] • 
In this case, the dynamics is described by four evolution 
equations for the interface height, electric field, and lipid 
densities of the two monolayers. The ability to track the 
lipid density while the interface is deforming provides 
valuable additional information. We find that during de- 
formation the density of the lipids on the top and bottom 
interface decreases in tandem, implying poration, see Fig- 
ure [I] Full details can be found in [25] . 




FIG. 4: (Color online) Increasing shape undulation in response 
to a perturbation in the transmembrane potential. The arrows 
show the fluid velocity field at t = 0. The color coding reflects the 
magnitude of lipid density where lighter shades indicate a decrease 
in the density (e.g the center region). Parameters as in Figure p] 

In conclusion, we have shown that the combined action 



of capacitive current through the membrane and mis- 
match of bulk conductivities destabilizes the membrane. 
This mechanism for creating shape undulations may be 
relevant for budding and vesiculation during electrofor- 
mation. Moreover, considering a membrane modeled as 
two coupled monolayers, we find that the lipid density de- 
creases in the regions of maximal deformation, which in- 
dicates poration in the context of the linearized model. A 
nonlinear investigation is needed to further clairify these 
results. 
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